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Spatially Resolving Substructures within the Massive Envelope 
around an Intermediate-mass Protostar: MMS 6/OMC-3 

Satoko Takahashi/ Kazuya Saigo,^ Paul T.P. Ho^'^ and Kengo Tomida,^'^ 

ABSTRACT 

With the Submillimeter Array (SMA), the brightest (sub) millimeter con- 
tinuum source in the OMC-2/3 region, MMS 6, has been observed in the 
850 fim continuum emission with approximately 10 times better angular resolu- 
tion than previous studies (~0.3"; ;^120 AU at Orion). The deconvolved size, the 
mass, and the column density of MMS 6-main are estimated to be 0". 32x0". 29 
(132 AUxl20 AU), 0.29 Mq, and 2.1x10^5 cm'^, respectively. The estimated 
extremely high mean number density, 1.5x10^° cm~'^, suggests that MMS 6-main 
is likely optically thick at 850 /xm. We compare our observational data with three 
theoretical core models: prestellar core, protostellar core + disk-like structure, 
and first adiabatic core. These comparisons clearly show that the observational 
data cannot be modeled as a simple prestellar core with a gas temperature of 
20 K. A self-luminous source is necessary to explain the observed flux density 
in the (sub) millimeter wavelengths. Our recent detection of a very compact and 
energetic outflow in the CO (3-2) and HCN (4-3) hues, supports the presence 
of a protostar. We suggest that MMS 6 is one of the first cases of an interme- 
diate mass protostellar core at an extremely young stage. In addition to the 
MMS 6-main peak, we have also spatially resolved a number of spiky structures 
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and sub-clumps, distributed over the central 1000 AU. The masses of these sub- 
clumps are estimated to be 0.066-0.073 Mq, which are on the order of brown 
dwarf masses. Higher angular resolution and higher sensitivity observations with 
ALMA and EVLA will reveal the origin and nature of these structures such as 
whether they are originated from fragmentations, spiral arms, or inhomogeneity 
within the disk- like structures/envelope 

Subject headings: stars: individual (0MC3-MMS 6) — stars: formation — ISM: 
evolution — ISM: clouds 



1. INTRODUCTION 

The study of the earliest phase of star formation, especially around the prestellar core 
and Class phases, is essential to understand star formation. This period includes the 
important phases such as the prestellar core phase, first adiabatic core phase, and the subse- 
quent formation of the protostar (Larson 1969; Masunaga & Inutsuka 2000, Bate 1998; 2010; 
Machida et al. 2006; Saigo et al. 2008; Schoneke & Tscharnuter 2011). Studying prestellar 
cores would provide clues for understanding the initial conditions of star formation. First 
adiabatic cores are able to launch outflows (Tomisaka 1998; 2002; Machida et al. 2006; 
Tomida et al. 2010a), and may be related to the formation of binary systems (Bonnell 1994; 
Matsumoto et al. 2003; Commercon et al. 2010). In the protostellar phase (i.e., after 2nd 
collapse), substantial mass accretes onto the central star. Hence this phase would be a key 
to determine the final stellar mass. However, candidates of this youngest stage of protostars 
are still very limited in number, and theoretical models are not well constrained as yet. 

The Orion Molecular Cloud-2/3 is located at the northern part of the Orion giant 
molecular cloud A {d^AlA pc; e.g., Menten et al. 2007), and is one of the nearest active star 
forming regions (OMC-2/3; e.g., Aso et al. 2000; Chini et al. 1997; Lis et al. 1998; Johnstone 
& Bally 1999; Nielbock et al. 2003; Takahashi et al. 2006; 2008a,b 2009; Peterson & Megeath 
2008). Previous studies identified approximately 50 individual star- forming sites. Detected 
sources are deeply embedded within the molecular filaments and most of these detected 
sources are identified as prestellar cores or protostellar cores (e.g., Chini et al. 1997; Lis et 
al. 1998; Johnstone & Bally 1999; Nielbock et al. 2003; Takahashi et al. 2006; 2008a, b; 
Peterson & Megeath 2008). 

With intcrfcromctric continuum observations, Matthews et al. (2005) and Takahashi 
et al. (2009) found a bright continuum peak (hereafter MMS 6-main). which has at least a 
factor of five larger fiux density at (sub) millimeter wavelengths as compared to all the other 
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OMC-2/3 sources (Takahashi et al. 2009). However, no star formation signatures such as 
a clear CO outflow, a radio jet, or an infrared source (<8 /um), have been detected at the 
source center in the previous studies (Matthews et ah 2005; Takahashi et al. 2009). The core 
density profile suggests a power-law structure with an index of -2.0. This power-law index 
suggests that the core still retains the density structure of the prestellar phase (Takahashi et 
al. 2009). The bolometric luminosity and core mass of MMS 6, as derived from the previous 
single-dish millimeter to far- infrared data, are <60 Lq and 36 M0, respectively (Chini et 
al. 1997). These parameters are approximately one order of magnitude higher than the 
values derived in typical low-mass star-forming regions such as the Taurus molecular cloud 
(c.f., Myers & Benson 1983). Therefore, MMS 6 would be one of the first cases for the 
detailed study of the intcrmcdiatc-mass star forming core at the earliest evolutionary stage 
(c.f., Beltran et al. 2008; Fuente et al. 2007; Takahashi et al. 2009; Enoch et al. 2009). 

In order to study the physical properties and the evolutionary stage of MMS 6-main 
more accurately, we have obtained the highest angular resolution image possible with the 
SMA, utilizing all the configurations at 850 //.m. Our observations have achieved approx- 
imately 10 times better angular resolution in terms of beam area, as compared to that of 
the previous SMA observations by Takahashi et al. (2009). We compare our results with 
the theoretically predicted star forming core models such as the prestellar core model, the 
protostellar core plus central disk model, and the first adiabatic core model. We make our 
comparisons in the visibility- amplitude domain and also in terms of the spectral energy 
distributions (SEDs). 

We describe the observations in Section 2. Observational results including sub-arcsecond 
images, physical properties, and core internal structures, will be presented in Section 3. 
Comparisons between the observed core density structures and the theoretical models, as 
well as the evolutionary stage of MMS 6-main, will be discussed in Section 4. Section 5 
summarizes the paper. 

2. OBSERVATIONS AND DATA REDUCTION 

The observations have been performed with the SMA (Ho et al. 2004), utilizing the 
sub-compact, compact, extended, and very extended configurations. The 850 /im continuum 
data were taken with seven antennas in the extended configuration (EX; projected baseline 
ranging between 40 and 240 kA), and the very extended configuration (VEX; projected 
baseline ranging between 90 and 580 kA) on February 9 and September 2, 2010, respectively. 
The typical system noise temperatures in DSB mode were between 200-350 K for both 
the EX and VEX configurations at the observed elevations. Both the LSB and USB data 
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were obtained simultaneously with the 90° phase switching technique by the digital spectral 
correlator, which has a bandwidth of 4 GHz in each sideband. After subtracting the channels 
where molecular lines might be expected within the band, the LSB and USB continuum data 
separated by 10 GHz RF were combined to improve the sensitivity. The effective bandwidth 
for the continuum emission is approximately 7.5 GHz. The phase and amplitude calibrator, 
0423-013 (1.5 Jy for EX and 2.9 Jy for VEX), was observed every 18 minutes. Observations 
of Callisto (for EX) and Vesta (for VEX) provided the absolute scale for the flux density 
calibrations. The overall flux uncertainty was estimated to be ~20%. The passband across 
the bandwidth was determined from observations of 3C 454.3 with a ~30 minute integration. 
The EX and VEX configuration data were combined with the similar wavelength continuum 
data taken with the SMA sub-compact and compact configurations. Detailed observing 
settings of the sub-compact and compact configurations are described in Takahashi et al. 
(2009; 2012). 

The raw data were calibrated using MIR, originally developed for the Owens Valley 
Radio Observatory (Scoville et al. 1993) and adopted for the SMA. After the calibration, 
the data from all the configurations were combined using the AIPS task "DBCON", and 
final CLEANed images were made using the DIFMAP task. In order to reduce residual 
phase errors and improve the dynamic range of the images, self-calibrations have been ap- 
plied. This improved the maximum dynamic range by a factor of 4. The resulting synthe- 
sized beam sizes of the combined data set were 0". 35x0". 30 (or 145 AU x 124 AU at the 
adopted distance of 414 pc) with a position angle of 32° for the uniform weighting image, 
and l".9xl".7 (or 787 AU x 704 AU) with a position angle of 76° for the Gaussian tapered 
image (FWHM=60 kA). The achieved rms noise levels were 2.7 mJy beam~^ for the uni- 
form weighting image and 12 mJy beam~^ for the Gaussian tapered image. Note that the 
self-calibrated images lose their absolute positions. The positional offsets between the non- 
self-calibrated image and the self-calibrated images (i.e.. Figure la, b) are (Tsl~0".048, and 
0".12 for the uniformly weighted and the tapered images. 

In addition, the positional accuracy is mostly determined by the phase drift due to the 
baseline error, ((Jbs~0.1A~0".1) of the non-self-calibrated image B Finally the positional 
accuracies of the final CLEANed images are estimated to be c7totai~\/o"BS^ + crsL^~0".ll for 
the uniform weighted image and ?a0".16 with the Gaussian tapered image. 



^The baseline error is also measured using the primary calibrator (0423-013) against the secondary cali- 
brator (0530-1-135). Both calibrators are located ~18° from MMS 6-main, and the separation between two 
calibrators is ^-^23°. The positional shift from the phase center of the secondary calibrator was measured to 
be ~0".l, which is consistent with the expected baseline error of the SMA, O.IA. 
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3. RESULTS 

3.1. Subarc-second 850 fim Continuum Image 

Figure 1 a and b sliow tlie 850 fim continuum images, which are predominantly thermal 
dust emission (Takahashi et al. 2009), obtained with the SMA. As denoted by crosses in 
Figure la, five infrared sources (i.e., lRSl-5; protostellar and young stellar objects) are 
located around MMS 6-main, but none of them are associated with the peak position of the 
850 fim continuum emission of MMS 6-main. Note that the positional shift between MMS 6- 
main and the closest infrared source, IRS 3, is 2". As discussed in Section 3.1.2 in Takahashi 
et al. (2009), the accuracy of the positional alignment between the infrared image and 850 
/xm image is on the order of 0".3, which was estimated from a common compact continuum 
source obtained in the same region and during the similar observational period. In addition, 
the positional accuracy of MMS 6-main obtained with the SMA 850 fim (uniformly weighted 
image) is 0".ll. These comparison show that the positional shift between MMS 6-main and 
IRS 3 is most likely real. 

In order to derive the peak fiux, fiux density, and source size, a 2D-Gaussian fit with 
two components, was applied to the MMS 6-main map. The derived properties for a com- 
pact and an extended component, are reported in Table 1. The residuals after the two- 
component fit are the 20% level at most as compared to the peak fiux. The non-uniform 
gas distribution, with substructures as described in the later part of this Section as well 
as non-Gaussian shaped envelope structures, are the causes of the residuals. The total 
fiux and the peak intensity of the compact component are measured as 1.1 Jy and 0.56 
Jy beam"^, respectively. The brightness temperature is estimated to be 52 K from the 
peak intensity within the beam size. The gas temperature is equal to the brightness tem- 
perature when the gas is optically thick (T^ = Tgas for r>l), while the gas temperature is 
larger than 52 K when the gas is optically thin (T^ = rTg^s for r << 1). Assuming the gas 
temperature is 52 K, r=l, the gas-to-dust ratio of 100, the column density is estimated to 
be A^H2 = (2.0±0.4)xl025 cm'^. This gives us the source mass of M=0.28±0.06 M©, and the 
mean number density of n//2 = 1.4±0.3xl0^° cm~'^, with the assumption of a spherically sym- 
metric structure §. For these estimates, the dust emissivity index of /9=0.93 (Takahashi et 
al. 2009) and the dust absorption coefficient of ka=0.027 cm^ g~^(A/400 /im)~^ are assumed 
(Keene et al. 1982). 

The total SMA fiux density of the compact continuum component associated with MMS 



•^Fitting errors of the deconvolved size and flux as well as the 20% absolute flux uncertainty described in 
Section 2 were considered in order to estimate the error of physical parameters. 
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6-main corresponds to 17 % of the flux density measured within the JCMT/SCUBA 13" 
beam at 850 fim (6.6 Jy; Johnstone & Bally 1999). Note that the ratio in beam surface area 
between the JCMT and the SMA is ~1900. A large amount of flux at the submillimeter 
wavelength, is concentrated within the central 100 AU size scale. 

In addition to the MMS 6-main component, a number of substructures extending from 
the MMS 6-main peak, have been spatially resolved in our continuum image (denoted as 
dashed blue arrows in Figure Ic and Figure 2). Here, we identify "substructures" as those 
distinct features whose peak fluxes have a signal-to-noise ratio of more than 25o", and whose 
contrasts with respect to their fainter surroundings arc more than lOo". We argue that these 
substructures, at the level of 10 times fainter than the MMS 6-main peak, are not artifacts 
due to sampling or image reconstruction. This is verified by calibrating and imaging our 
secondary cahbrator, 0530-1-135, in exactly the same manner as for MMS 6. In this image, a 
dynamic range of at least 50 is achieved, with no hint of any substructure for the secondary 
calibrator, down to the theoretical noise level. We have also compared the maps made with 
uniform or natural weighting, and with or without self-calibration. The same structures are 
seen with good correspondence between the different maps. Hence we conclude that the 
detected substructures are robust. 

Among the detected substructures, four are particularly well resolved. We mark them 
with orange circles in Figure Ic and Figure 2. Their contrast with the surrounding structures 
at the same radii is much greater than 10 sigma. We call these four, the "sub-clumps". All 
these sub-clumps arc located at 0".9-l".7 (370-700 AU) from the MMS 6-main peak position. 
The sizes of the sub-clumps are comparable to the SMA synthesized beam. 

The peak fluxes of the sub-clumps are 70-80 mJy beam"^. Assuming that the dust 
emission is optically thin within these sub-structures, and that the temperature distribution 
of the dust continuum emission is uniform with Td=20 K, the masses of the sub-clumps 
are estimated to be M=0.066 - 0.073 M0. Such masses are on the order of the values for 
brown dwarfs, (0.011<M*<0.084 M©). Assuming that the sub-clumps have a similar size 
as the synthesized beam and a spherical geometry, the number densities are estimated to 
be nij2f=:i3xl0^ cm~^. Figure 2 shows the flux distribution along the azimuthal direction as 
a function of radius from the MMS 6-main peak position. The flux fluctuations along the 
azimuthal directions, are exactly what we had noted in Figure Ic. Separations between these 
sub-clumps/spikes range between 0.6-1.4" (250-580 AU; as denoted by arrows in Figure 2). 
The origin of these sub-structures is discussed in Section 4.3. 

In addition, extended emission at the 5a level, is detected to the southeast of MMS 6-main. 
This elongation is consistent with the large scale filamentary structure (e.g., Chini et al. 
1997), and may reflect the fragmentation processes within the filaments. 
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3.2. Visibility-Amplitude Plot 

The high flux density concentration is attributed to not only the density distribution, 
but also to the temperature distribution. In order to deflne these distributions closer to 
the center, we further resolve the core internal structure as first elucidated by Takahashi 
et al. (2009). We use the newly obtained SMA EX and VEX configuration data to study 
approximately one order of magnitude smaller size scales {R ~60 AU) as compared to the 
previous analysis. Core internal structures are discussed in the Fourier domain analysis 
by utilizing the visibility-amplitude plot as presented in Figure 3. We directly analyze the 
interferometric data, bypassing the non-linear image deconvolution process and the problem 
of the missing fiux. 

Note that the binned amplitude-visibility plot is processed with the vector amplitude 
average. Substructures shifted from the phase center of the visibility plot will be atten- 
uated due to the coherence loss. This effect suppresses the ffux contamination from the 
substructures. Asymmetric structures within the core are presented as a variation of the 
visibility-amplitudes as a function of the projected wf-distance. 

We assume that the density and temperature follow simple radial power laws as p(r)ocr~^ 
and T{r)c<.r~'^, respectively. Then the power law of the visibility amplitude plot, 7, have a 
relation of'~f = p + q — 3 (e.g., Harvey et al. 2003a; Takahashi et al. 2009). Figure 3 shows 
the visibility amplitude plot as a function of the uv distance. Our observational results 
suggest two different power-law indices, with the break at 60 kA {D =1300 AU). The outer 
part and inner part of core are fitted by 7 =-0.64 and -1.4, respectively. In the isothermal 
case (q—O), the density power- law index of the outer and inner part of core are derived to be 
2.4 and 1.6, respectively. Assuming the central heating source (g=0.41 derived by Takahashi 
et al. 2009; optically thin case), the density power-law index of the outer and inner part of 
core are derived to be 2.0 and 1.2, respectively. These results suggest that the density profile 
is shallow at the inner part of core (D<1300 AU). 



4. DISCUSSION 

4.1. Model Compcirisons 

In order to further discuss the internal structure, the observed results are compared with 
core models in the different evolutionary stages based on the visibility-amplitude (Figure 3) 
and the spectral energy distribution (Figure 4): (i) Prestellar core, (ii) First adiabatic core, 
and (iii) Protostellar core. For all models, the gas temperature of 20 K, which is the typical 
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gas temperature derived from the northern part of OMC clumps, is adopted (Cesaroni & 
Wilson 1994). Moreover, an effective sound speed of 0.42 km s~^, which is derived from the 
H^^CO"'" observations (Takahashi et al. 2009), was included in order to take account of the 
observed non-thermal motions. The radiative transfer calculations are performed to produce 
the visibility-amplitude plot and SED for model cores. For the calculations, composite 
aggregates dust model (Semenov et al. 2003) was adopted. 

(i) Prestellar core: Larson-Pcnston model (ccfr=0.42 km s~^ and t = 0; the end of the 
prestellar phase) is adopted (denoted by a blue curve in Figure 3). Note that the Larson- 
Penston model is dependent only on the sound speed. The measured sound speed of 0.4 
km s~^ by Takahashi et al. (2009) was adopted. The simple prestellar core model shows 
approximately one order of magnitude smaller visibility amplitude distribution as compared 
with the observational result. The results unlikely show that MMS 6- main is the prestellar 
core. Recent discovery of an extremely compact molecular outflow associated with MMS 
6-main by Takahashi & Ho (2012) supports to discard the prestellar scenario. 

(ii) The first adiabatic core: This model, denoted by a green curve in Figure 3, is 
calculated using the three dimensional radiative hydrodynamic numerical simulations (e.g., 
Tomida et al. 2010b; Saigo et al. 2012). In order to investigate the formation of the interme- 
diate mass star, we adopt the initial condition of the rotating critical Bonnor-Ebert sphere 
with the effective sound speed of Ccs =0.42 km s~^, which was measured from the H^'^CO"'" 
(1-0) observations (Takahashi et al. 2009). The central density of n//2=2.6xl0^ cm^'\ which 
is a typical density for nearby molecular clouds as well as the value usually used for the first 
adiabatic core simulations (e.g., Saigo et al.2008; 2011), was adopted as an initial density for 
the parental molecular cloud. In this calculation, we assume the gas temperature of 20 K. 
The age of the plotted flrst adiabatic core is f=:i2000 yr. The plotted model is the best model 
to reproduce observed flux density among their evolution. Resulting visibility amplitude dis- 
tribution are denoted by the green curve in Figure 3. The shape of the visibility distribution 
is not consistent with the observed data. Note that younger phase of the simulated first 
adiabatic core models show lower flux values than presented model due to less significant 
gravitational potential of the first core, while later phase of the first adiabatic core models 
does not show much differences in terms of the total flux, but show even flatter visibility 
amplitude distributions for much of the baseline range compared to the presented model. 

(iii) Protostellar core: This model assumes the Larson-Penston core + central disk 
structure (denoted by a red curve in Figure 3). Parameters of the protostellar model are 
selected for reproducing the observed amplitude. For the disk-like structure, R = 300 AU, 
M — 1.5 M©; S(r)cxr~^, and T=30 K are assumed. Moreover, the temperature distribution 
in the envelope of the central protostar was assumed as T(r) = 500(r/10 AU)~^ K. The given 
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temperature gradient, T{r)(x.r ^, is considering a central heating source and an optically 
thick envelope (Hartmann 1998). 

The prestellar core model has a lower flux as compared with the other two models. This 
is due to the cool gas temperature (20 K) as expected in a prestellar core. The observed 
large flux at 850 //m is not consistent with the prestellar core model, but clearly suggests 
that a self-luminous source (i.e. higher temperature gas than 20 K) is required to explain 
the observed large flux at 850 //m. 

Note that the visibility amplitudes present a bump or excess between 200 kA and 400 kA 
as denoted by the blue arrows in Figure 3 as denoted by a blue arrow. However, the bump 
observed at Ri200 kA corresponds to the linear size of i?f«300 AU. We can produce such a 
bump in the protostellar core by adding a 30 K, 300 AU disk-like structure as we shown 
in protostellar model in Figure 3. The addition of such a disk is too simplistic, in view of 
the observed structures. The observed bump in the visibility amplitudes above a smooth 
extended core, is due to the clumpy structures including the MMS 6-main peak which we 
have detected in Figure Ic, and which all have a characteristic size scale of around 300 
AU. Furthermore, the visibility amplitude variations are a function of the azimuth angle 
(Figure 2). 

Figure 4 presents the comparisons between the observed MMS 6-main spectral energy 
distribution and their core models. Black open squares in Figure 4 show the observed SED of 
MMS 6-main measured with an angular resolution of 0".6-6."0. These values were retrieved 
from Table 5 in Takahashi et al. (2009). Spitzer/MIPS 24 //m and 70 //m upper limits are 
also shown in the SED. These upper limits are due to the contamination from another Class I 
source, IRS 3, located at 2" north of MMS 6-main. On the other hand, 850 /im continuum 
emission obtained with the SMA, 2.3 mm and 3.3 mm continuum data obtained with the 
Nobeyama Millimeter Array, and 7.3 mm continuum emission obtained with the Very Large 
Array, are from lower limits due to the missing flux. However, there are no contamination 
from nearby sources at these wavelengths because of the high-angular resolution observations. 

Observed millimeter/submillimeter continuum emission have higher flux densities as 
compared with the SED expected in the prestellar and first adiabatic core models denoted 
by the blue line and the green line in Figure 4. For example, approximately a factor of 
six higher flux density was observed at 850 //m as compared to these two models. The 
protostellar model as denoted by the red hue agrees well with the observed flux densities. 
SED comparisons clearly imply that the observed high flux densities are caused by a self- 
luminous source. 
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4.2. Evolutionary Stage of MMS 6-main 

In terms of the core, the sizes of MMS 6-main and the other protostellar sources in the 
OMC-3 region do not vary significantly (i.e., core sizes measured at more than Scr signal 
level at 850 /xm range between 900 AU and 2100 AU; Takahashi et al. 2012). This suggests 
that MMS 6-main may have a similar density structure in the outer most region. However, 
MMS 6-main has more than a factor of five larger millimeter/submillimeter flux densities as 
compared with other sources detected in OMC-2/3 with any aperture size between 10000 AU 
and 500 AU (e.g., Chini et al. 1997; Takahashi et al. 2009; 2012). 

The visibility amphtude plot and SED clearly imply the presence of a central heating 
source embedded within MMS 6-main. Hence, MMS 6-main is most likely in the protostel- 
lar core phase. The absence of a mid-infrared source, a large-scale molecular outflow, or 
an ionized jet, with current sensitivity and angular resolution, can be explained if the cen- 
tral driving source is extremely young (less luminous with high-extinction from surrounding 
circumstellar material), so that the outflow and jet are not yet substantial. Due to the 
high-extinction, the infrared emission is absorbed by the massive surrounding envelope and 
then reradiated in the submillimeter wavelengths. Takahashi & Ho (2012) have recently dis- 
covered an extremely compact molecular outflow (1000 AU) associated with MMS 6-main 
in CO (3-2) and HCN (4-3). This is the smallest molecular outflow associated with the 
intermediate-mass protostars and also even smaller than the outflows found in the proto- 
stellar cores in the Taurus star forming region (e.g., Hogerheijde et al. 1998). Nevertheless, 
the estimated outflow force (~10~^ M0 km s~^ yr~^) is similar as compared to the other 
large-scale molecular outflows associated with the protostellar cores in the OMC-2/3 re- 
gion. Furthermore, the observed bow-shock type of outflow velocity structures are similar 
to those observed in the Class type protostars. These results support our suggestion that 
MMS 6-main is in the protostellar phase, and in particular at a very early evolutionary 
phase. 



4.3. Substructure within the Massive Core 

Substructures were detected in our sub arcsecond resolution of 850 fim continuum image 
for the first time. There are two possible formation mechanisms to produce the substruc- 
tures and sub-clumps surrounding MMS 6-main: (i) gravitational instability within 
dynamically collapsing cores, or (ii) gravitational instability within rotationally 
supported disks. 
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(i) dynamically collapsing core case: Hanawa & Matsumoto (2000) suggest that 
collapsing cloud cores become more unstable if the effective sound speed (e.g., turbulent mo- 
tions) decreases during the collapse (see also Saigo et al. 2000). Non-thermal velocity width 
of 0.42 km s~^ (mach=2.5) are observed in the outer part of MMS 6 core (~1500 AU). The 
Jeans length is described as Afrag = \/7rc^/Gpo- Here, Cg is the sound speed, and the relation 
between the sound speed and the temperature of gas (T^) is described as Cg = ^JkTd/ ^imu- 
G, po(= l^i^HnH2)i 1^1 cind rriH are the gravitational constant, the mean density of the gas, 
and the hydrogen mean molecular weight of 2.33, and the hydrogen mass, respectively. Note 
that both the temperature and density vary slowly and in the same sense, so that the Jeans 
analysis is approximately correct as an estimate. Adopting the gas temperature of 20 K, 
and the typical number density of —3.0x10^ cm~^ (derived in Section 3.1), the Jeans 
length is estimated to be 360 AU (0."87). Note that the number density derived from the 
sub-clumps is used to represent the number density of the parental gas (i.e., the number 
density before the fragmentation). Hence, the estimated Jeans length is the lower-limit of 
the fragmentation scale, if fragmentation occurred much earlier. The derived Jeans length 
is similar to the separation between the sub-clumps (i.e., the separations of local peaks as 
measured in Figure 2), which is 0.6"-1.4" (250-580 AU). 

(ii) rotationally supported disk: Shu et al. (1990) show that if the rotationally sup- 
ported disk becomes sufficiently massive, Md/{M^. + Mdis)~24%, the rotationally-supported 
disk becomes gravitationally unstable and spiral arms are formed. In the MMS 6-main case, 
a 2.5 Mq disk, which corresponds to Md/{M^ + Mdisk)~45%0, might be suggested from the 
observations. In such a case, it is possible that the observed sub-structures could be due 
to instability within the purported disk. Fragmentation of self-gravitating massive disk-like 
structure has also been suggested by other studies (e.g., Gammie 2001; Kratter & Matzner 
2006; Kratter et al. 2008). We emphasize that currently, our observations only indicate 
substructures, without any definitive evidences of a disk-like structure. With current sensi- 
tivity, we also do not have any kinematic information from the envelope/disk-like structures. 
Sensitive molecular line observations are required to investigate the kinematics and to test 
the proposed scenario. 

In both cases, the physical processes follow the self-gravitational forces. Moreover, it 
is important to note that fragmentation within the disk will disturb the mass accretion 
processes onto the central star (e.g., Kratter et al. 2006; 2008). This process plays an 
important role in determining the final mass of the primary star. Our observations are 



''Here, stellar mass of 3 Mq was assumed as a maximum stellar mass at this stage. Mass of the disk-like 
structure was estimated from the total flux density detected at more than 5a level (4.7 Jy), after subtracting 
the flux contribution from MMS 6-main (2.0 Jy) with the assumption of Td=2Q K 
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the first case to present fragmented sub-structures within core/disk-hke structure around 

an embedded protostar with i?^500 AU. As we mentioned in Section 3.1, the mass of sub- 
clumps (0.066-0.073 Mq) have brown dwarf mass. Origin of the detected sub-clumps could 
be the eventual cluster members around the primary star or related to brown dwarf/planet 
formation. 



4.4. Future Prospects 

Limited data quality at the projected uv distances of longer than 400 kA, as well as 
the presumably optically thick continuum emission at 850 /xm, make it difficult to constrain 
the detailed spatial structures within Z}<100 AU scale. Even higher angular resolution 
observations with the relatively optically thin wavelengths (A>850 jim. ) will be necessary 
for defining the detailed internal density and temperature structures. Observations with 
the EVLA (e.g, K, Ka, and Q band) and ALMA (e.g., Band3-6) will be a key for these 
studies. These will also reveal detailed structures within the spikes and sub-clumps and 
observationally constrain the origin of these structures. Observations using the high density 
tracer or hot core tracer with the SMA, ALMA and EVLA are important for revealing gas 
physical conditions and dynamics of the inner envelope, as well as the circumstellar disk (i.e., 
rotationally supported disk), if present, at the earliest evolutionary phase (e.g., Zapata ct al. 
2010). Moreover, high angular resolution (<< 1") and high sensitivity EVLA observations 
are also crucial in order to search for the mass ejection phenomena associated with the central 
star such as faint ionized jet and/or H2O maser. 

5. Conclusion 

With the highest angular resolution possible with the SMA, an intermediate-mass pro- 
tosatar, MMS 6-main, has been spatially resolved. The observed flux density at the MMS 
6-main (Dail20 AU) corresponds to the brightness temperature of 52 K. The column density 
and the mass are estimated to be Nh^ =2.1x10^^ cm~^ and M=0.29 M©, respectively. De- 
tailed model comparisons clearly show that the observed visibility amplitudes cannot be ex- 
plained by the simple prestellar core model with a gas temperature of 20 K. A self-luminous 
source is necessary to explain the observed high- visibility amplitudes. Multi- wavelength 
comparisons using the SED as well as recent extremely compact molecular outflow detec- 
tion by Takahashi & Ho (2012) also support a protostellar type source. Several sub-clumps 
and spikes have been spatially resolved. These structures are well centered at MMS 6-main. 
Masses of the detected sub-clumps are estimated to be on the order of brown dwarf mass 
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(0.066-0.073 Mq). The origin of the density fluctuations can be explained by the gravita- 
tional instability within the massive core or a massive disk-like structure. These structures 
could form brown dwarfs or planets around the primary star. 
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Table 1. Parameteres of MMS 6-main 





RA 

[J2000] 


Dec 
[J2000] 


Deconvolved Size, P.A. 
[arcsec, degree] 


Peak Intensity 
[Jy beam~^] 


Flux Dei 

[Jy] 


Compact 
Extended 


05 35 23.42 
05 35 23.41 


-05 01 30.57 
-05 01 30.57 


0.318±0.004x0.295±0.004 (174±7°) 
1.788±0.001xl.448±0.001 (156±2°) 


0.556±0.003 
0.152±0.002 


1.061±0 
3.931±0 




r. . I . , . . I . . , . I . . , . I . . . , I . H M" --- L^^ 

5'35~34.5" 34.0' 23.5" ss.o" 23.5" s''36"'33,7'23.6' 3S,S' 23.4' E3.3' Si^S^S 3.7" 23.8' 23.5" 23.4" 23.3= 



Right .isfenslon (J^OOO) 

Fig. 1. — The 850 fim continuum images taken with the SMA subcompact, compact, ex- 
tended, and very extended configurations, (a): Image made with uv tapering of FWHM=60 
kA. The contour levels are ±5a, 10a, 15a, 20a, 25a, 30a, 40a, 60a, 80a, 120a, 160a, 200a, 
and 240a (la=12 mJy beam~^). (b): Image made with uniform weighting. The contour 
levels are ±5a, 10a, 15a, 20a, 25a, 30a,40a, 60a, 80a, 120a, 160a, 200a, and 240a (la=2.7 
mJy beam~^). (c): Same as figure (b), overlaid with notes for detected substructures as 
denoted by blue dashed arrows and the sub-clumps as denoted by orange dashed circles. 
Denoted numbers for each substructure are corresponding to same numbers as in Figure 2. 
Negative contours are presented by the dashed lines. The crosses show the positions of the 
infrared sources identified in Takahashi et al. (2009). The size of the crosses corresponds to 
the positional accuracies of the Spitzer images in Takahashi et al. (2009). Filled ellipses at 
the bottom-left corners show the synthesized beams of each image. These figures were made 
using the KVIS in KARMA (Gooch 1995). 
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Fig. 2. — The flux distribution as a function of azimuthal direction and radius centered at 
the peak position of MMS 6-main. 1 pixel corresponds to ff' .03. The black lines in the 
vertical axis show 1" and 1." 5 from the center. Positions of spiky substructures are denoted 
by blue dashed allows. Sub-clumps denoted by dashed orange circles. Denoted numbers of 
each substructure are corresponding to same numbers denoted in Figure Ic. East direction is 
defined as the angular offset as (F and is counted toward counter clock direction as a positive 
offset. Figure was made using the KPOLAR in KARMA (Gooch 1995). 
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UVdistance (k lambda) 

Fig. 3. — Visibility amplitude plot binned by 4 cls a function of the projected uv distance 
produced from the 850 ^m continuum data. Black open circles show the observational results 
obtained from the SMA 850 fim continuum image. The statistical errors (la) of the 850 fim 
data are denoted by vertical error bars. For the comparisons, prestellar core model (Larson- 
Penston core) with 20 K gas is denoted by blue curve. First adiabatic core model is denoted 
by green curve. Protostellar core model (Larson-Penston core with disk-like structure) is 
denoted by red curve. The core models are calculated with the inclination angle of 3(F . 
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Fig. 4. — Spectral energy distributions measured in the millimeter, submillimeter, and mid- 
infrared wavelengths are denoted by open squares. These data were retrieved from Takahashi 
et al. (2009). The bars within the open squares show the error bar of the absolute flux 
measurements (i.e., ±20% for 850 nm, ±15% for 2.3 mm and 3.3 mm, and a few percent for 
1.3 mm). Model SEDs were overlaid with the observed data: prestellar core model (Larson- 
Penston core) with 20 K gas is denoted by blue curve. First adiabatic core model is denoted 
by green curve. Protostellar core model (Larson-Penston core with disk-like structure) is 
denoted by red curve. The core models are calculated with the inclination angle of 30° . 



